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Abstract 



A search for microscopic black holes and string balls is presented, based on a data 
sample of pp collisions at = 8 TeV recorded by the CMS experiment at the Large 
Hadron Collider and corresponding to an integrated luminosity of 12 fb -1 . No ex- 
cess of events with energetic multiparticle final states, typical of black hole produc- 
tion or of similar new physics processes, is observed. Given the agreement of the 
observations with the expected standard model background, which is dominated by 
QCD multijet production, 95% confidence limits are set on the production of semi- 
classical or quantum black holes, or of string balls, corresponding to the exclusions of 
masses below 4.3 to 6.2 TeV, depending on model assumptions. In addition, model- 
independent limits are set on new physics processes resulting in energetic multipar- 
ticle final states. 



Submitted to the Journal of High Energy Physics 



© 2013 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 



*See AppendixJXjfor the list of collaboration members 



1 



1 Introduction 

Theoretical models with low-scale quantum gravity aim to account for the origin of the large 
difference between the electroweak scale (~0.1TeV) and the Planck scale (Mpj ~ 10 16 TeV), 
known as the hierarchy problem of the standard model (SM) of particle physics. One of the 
predictions of such scenarios is the possibility of producing microscopic black holes or their 
quantum precursors in proton-proton collisions at the CERN Large Hadron Collider (LHC) (H 

The basis of this analysis is the theoretical model proposed by Arkani-Hamed, Dimopoulos, 
and Dvali (ADD) (SHU. This model attempts to solve the hierarchy problem by introducing 
n large, flat, extra spatial dimensions, compactified on an n-dimensional torus or a sphere. 
By opening the multidimensional space only to the gravitational interaction, the fundamental 
Planck scale in 4 + n dimensions, Md, is lowered to the electroweak symmetry breaking scale, 
such that M D +2 o< Mp[R~" where R is the radius of the extra dimensions. The reduction in 
Md is accomplished without affecting tight constraints coming from precision measurements 
of properties of other types of fundamental interactions. The enhanced gravity in multidimen- 
sional space allows the formation of microscopic black holes. Production of black holes is also 
possible in the Randall-Sundrum (RS) model [5 - 7] and in models with unparticles DSHH. In the 
former case, the hierarchy problem is addressed by adding a single compact spatial dimension 
with radius comparable to the Planck length and with the metric of the 5-dimensional space- 
time being exponentially "warped" in the direction of this extra dimension (the anti-deSitter 
space-time). 

We present an extension of previous searches for microscopic semiclassical black holes (TJ |2|, 
quantum black holes ||5j[10|, and string balls [11] conducted by the Compact Muon Solenoid 
(CMS) Collaboration at the CERN LHC fl"2l|T3| . The analysis utilizes a data sample correspond- 
ing to an integrated luminosity of 12.1 ± 0.5 fb _1 , collected by the CMS detector in pp collisions 
at a centre-of-mass energy of 8 TeV. These data were recorded in the 2012 running period of the 
LHC. The 14% increase in the energy of the machine relative to the 2011 dataset would result in 
a larger production cross section for black holes and other new physics with energetic multipar- 
ticle final states. This allows the current search to penetrate a previously unexplored regime. 
Searches for black holes have also been performed by the ATLAS Collaboration iMllTBll . 

A characteristic signature of evaporating semiclassical black holes or string balls is a large num- 
ber of energetic final-state particles of various types, while quantum black holes typically decay 
into a pair of energetic jets. Further details of the analysis method and the underlying mod- 
els can be found in earlier publications |[T2l [13| . The results are presented in terms of a set 
of benchmark scenarios and are also interpreted in terms of model-independent limits on the 
production cross section for a new physics phenomenon multiplied by the branching fraction 
of its decay into a multiparticle final state. 

2 The CMS detector 

A detailed description of the CMS detector can be found elsewhere |16 |[ . The central fea- 
ture of CMS is a 3.8 T superconducting solenoid of 6 m internal diameter that encloses a sil- 
icon pixel and strip tracker, a lead-tungstate crystal electromagnetic calorimeter (ECAL), and a 
brass /scintillator hadron calorimeter (HCAL). Muons are detected in gas-ionisation detectors 
embedded in the steel flux return yoke of the magnet. 

The ECAL is a finely segmented calorimeter that uses crystals situated in a barrel region ( | tj \ < 
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1.48) and two endcaps that extend to \rj\ = 3.0. Here pseudorapidity rj is defined as —In [tan (6/2)], 
where 9 is the polar angle measured from the geometrical centre of the detector with respect 
to the anticlockwise proton beam. The transverse dimensions of the lead-tungstate crystals are 
Arj x Acp = 0.0174 x 0.0174, where cp is the azimuthal angle in radians. The HCAL consists of 
interleaved brass plates and scintillator sheets that extend to \t]\ = 3.0. The granularity of the 
HCAL towers is Arj x Acp = 0.087 x 0.087. Extensive forward calorimetry complements the 
coverage provided by the barrel and endcap detectors. 

The CMS trigger system is composed of two levels that are used to select potentially interesting 
events. The first level (LI) trigger ensures negligible dead time and is responsible for reducing 
the event rate to 100 kHz using the information from calorimeters and muon detectors. After 
LI triggering, the data are passed to the software-based high-level trigger, which decreases the 
event rate to several hundred Hz for further storage. 

The data used for this analysis are collected with a set of triggers based on the scalar sum of 
the transverse energies (Hj) of the calorimeter jets found by the trigger. The thresholds for the 
Hj triggers increased from 200 to 750 GeV, depending on the data-taking period, to cope with 
the increasing instantaneous luminosity of the LHC. For the earlier part of the data taking, 
we additionally utilized Hj triggers that use jets reconstructed using the particle-flow (PF) 
technique |fT7) , which are corrected for the calorimeter response to calculate the Hj variable. 
The trigger is measured to be fully efficient for jet-enriched collision events with Hj above 
1 TeV. 

3 Event reconstruction and Monte Carlo samples 

The PF technique is used offline to reconstruct and identify charged and neutral particles using 
information from all the subdetectors. Jets are reconstructed by clustering the PF candidates 
using an infrared-safe ardi-kj algorithm with a distance parameter of 0.5 ||T8Tl20l . Additional 
selection criteria are applied to jets to remove noise and non-collision background [21J. The PF 
jets are required to have transverse momentum pj > 50 GeV and to lie within \rj\ < 2.6. The 
jet energy response is further corrected using simulated events, as well as dijet and photon+jet 
collision events pi) . 

Muons are reconstructed using the PF algorithm by matching the tracks in the silicon detector 
to segments in the muon chambers. Muons with \t]\ < 2.1 and pj > 50 GeV are selected. 
Furthermore, they are required to have an impact parameter less than 0.2 cm to suppress the 
cosmic ray muon background. In addition, the scalar sum of charged and neutral particle 
transverse energies, calculated in a cone of AR = [Acp) 2 + {Arj) 2 = 0.3 around the muon 
direction, should not exceed 20% of the muon transverse momentum. The minimum separation 
between any two objects (jet, lepton, or photon) is required to be AR > 0.3. 

Electrons and photons depositing energy in the ECAL are identified via clustering algorithms, 
taking into account the expected cluster shapes. Electron reconstruction uses the PF algo- 
rithm ffl7ll22l and requires a silicon tracker trajectory to match an energy cluster found in the 
ECAL. Photon reconstruction uses ECAL clusters and requires no matching hits in the pixel 
tracker and an ECAL deposit with a shape consistent with that expected for a photon. Both 
objects are required to have pj > 50 GeV and to lie within the fiducial region of the barrel 
(\f]\ < 1.44) or endcap (1.56 < \tj\ < 2.4). The barrel-endcap transition region is excluded 
because the reconstruction of electrons and photons in this region is not optimal; energetic 
electrons and photons in this region nevertheless contribute to the reconstruction of jets. The 
separation AR between the electron candidate and any muon candidate that has more than 10 
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hits in the inner tracker is required to be greater than 0.1. We also require the scalar sum of 
charged and neutral particle transverse energies, calculated in a cone of AR = 0.3 around the 
electron direction, not to exceed 20% of the electron transverse momentum. The ratio of HCAL 
to ECAL energy deposits is required to be less than 5% for photon candidates. Photons must 
be isolated in the tracker, ECAL, and HCAL. The scalar sums of transverse energy (momenta 
in the case of the tracker) are calculated in a cone of AR = 0.4 around the candidate photon di- 
rection. These sums should not exceed 2.0, 4.2, and 2.2 GeV for the tracker, ECAL, and HCAL, 
respectively. 

The missing transverse energy (E™ 1SS ) is defined as the absolute value of the vector sum of 
transverse momenta of all the PF objects reconstructed in an event. The E™ 1SS measurement is 
corrected to account for the jet energy scale calibration ||23| . 

Simulated samples of semiclassical black hole events are produced using the parton-level Black- 
MAX v2.01 0241 and Charybdis vl.0.3 |25j[26j Monte Carlo event generators. Various models 
are simulated, including black holes that are produced non-rotating or rotating; those with or 
without mass and angular momentum loss at the time of formation; and those with or without 
a stable or boiling remnant. In addition, the modified BLACKMAX generator settings |27| are 
used to simulate the production of string balls. Detailed descriptions of each model can be 
found in Refs. HUSH. All these signal samples are generated using the MSTW20081o68cl (281 
parton distribution functions (PDF). Samples of quantum black holes are generated with the 
QBH vl.03 matrix element generator [29 1 using the CTEQ6L PDF set II30I . The parton-level 
events produced by these generators are then used as input to the pythia v6.426 |3T| parton 
showering simulation and a fast parametric simulation of the CMS detector [32, 33 j. The fast 
simulation was validated with the full detector simulation, using the Geant4 |34| framework 
for several benchmark points. 

The small backgrounds from 7 + jets, W + jets, and tt production are estimated from Monte 
Carlo simulations using the MadGraph v5 |35|1 matrix element event generator interfaced 
with the PYTHIA parton showering simulation, followed by the full detector simulation using 
Geant4. These background samples are generated using the CTEQ6L PDF set. 

4 Analysis method 

The analysis strategy is identical to the one used in the previous analysis at a/s = 7TeV and 
is described in detail in Ref. |[T3| . The search for black holes and string balls is performed 
using events with at least two jets and any number of photons and leptons. There is no explicit 
requirement of missing transverse energy in the event. 

The search for black holes is based on a search for a deviation from the SM background pre- 
dictions in the Sj spectra observed in data. The Sj variable is defined as the scalar sum of 
transverse momenta of all the final-state objects in the event with pj in excess of 50 GeV. If the 
E™ 1SS in the event exceeds 50 GeV, its value is also added to the Sj variable. We then determine 
the multiplicity N of the objects in the final state by counting all the objects in the events (ex- 
cluding E™ 1SS ) that enter the calculation of St |l2l[T3| . We analyse the data for various inclusive 
multiplicity bins, from N > 2 to N > 10, and look for deviations from the SM background 
predictions in each of these bins. 

The SM background is completely dominated by QCD multijet production and is estimated 
directly from data using a method based on Sj multiplicity invariance ffl2l[l3| . All other back- 
grounds are negligibly small in the Sj range used in this analysis, as shown in figure [TJ The 



4 
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multijet background estimation method is based on the empirical observation that the shape of 
the Sj spectrum is approximately independent of N, so the shapes of the St spectrum for any 
number of objects can be estimated using a fit to the dijet data (N = 2). The dijet sample has 
been previously studied in a dedicated analysis II36II and is known not to exhibit any signal-like 
features in the range of 1.8 < Sj < 2.8 TeV, which is used to obtain the background shape. The 
central value of the background shape and its uncertainty are determined from the fit to several 
semi-empirical template functions IT131 , by taking the best fit function as the central value and 
the envelope of the alternative fits as the measure of systematic uncertainty in the background 
shape. The background shape is parameterized with the function Pq(1 + x) Pl / x F2+Pilo ^ x \ and 
the uncertainty envelope is defined with two additional functions, Pq/(Pi + P^x. + x 2 ) P3 and 
Po/ (-Pi + x) Pl . Here, Pj are the fit parameters and x = Sj. Results of the fit can be seen in 
figure [T] 

The scaling of the background to higher multiplicities is performed by normalising the back- 
ground shape to data in each inclusive multiplicity bin in the control range (1.9 < Sj < 2.3 TeV), 
where any significant signal contribution has been already ruled out by earlier analyses |12l[T3| . 
The lower boundary of the control region is chosen to be substantially above the trigger and 
multiparticle (N x 50 GeV) turn-on regions. 

The Sj distributions for data, for predicted background, and for three semiclassical and one 
quantum black hole signal benchmarks are shown in figures [l]-[3] for a number of exclusive 
and inclusive multiplicities. We do not plot the quantum black hole signal St distributions for 
inclusive multiplicities of four or more, as the search for quantum black holes is not sensitive in 
higher inclusive multiplicity bins. No statistically significant excess of events over the expected 
background is observed in any of these spectra. 
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Figure 1: Distribution of the total transverse energy, St, for events with multiplicity: (Left) 
N = 2 and (right) N > 2 objects (photons, electrons, muons, or jets) in the final state. Ob- 
served data are depicted as points with error bars; the solid line with a shaded band is the 
multijet background prediction and its systematic uncertainty. Coloured histograms represent 
the 7 + jets (orange), W + jets (red), and tt (green) backgrounds. Also shown are the expected 
semiclassical black hole signals for three parameter sets of the BlackMax non-rotating semi- 
classical black hole model, as well as a quantum black hole model. Here, Mj^f is the minimum 
black hole mass, Mqbh is the minimum quantum black hole mass, Md is the multidimensional 
Planck scale, and n is the number of extra dimensions. 
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Figure 2: Distribution of the total transverse energy, St, for events with multiplicity: (Top left) 
N > 3, (top right) N > 4, (bottom left) N > 5, and (bottom right) N > 6 objects (photons, 
electrons, muons, or jets) in the final state. Observed data are depicted as points with error 
bars; the solid line with a shaded band is the multijet background prediction and its systematic 
uncertainty. Also shown are the expected semiclassical black hole signals for three parameter 
sets of the BlackMax nonrotating black hole model, as well as a quantum black hole signal 
of the QBH model. Here, Mj^f is the minimum black hole mass, is the minimum quan- 

tum black hole mass, Md is the multidimensional Planck scale, and n is the number of extra 
dimensions. 
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Figure 3: Distribution of the total transverse energy, Sj, for events with multiplicity: (Top left) 
N > 7, (top right) N > 8, (bottom left) N > 9, and (bottom right) N > 10 objects (photons, 
electrons, muons, or jets) in the final state. Observed data are depicted as points with error 
bars; the solid line with a shaded band is the multijet background prediction and its systematic 
uncertainty. Also shown are the expected semiclassical black hole signals for three parameter 
sets of the BlackMax nonrotating black hole model. Here, Mj^f is the minimum black hole 
mass, Md is the multidimensional Planck scale, and n is the number of extra dimensions. 
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5 Results 

In the absence of an excess of data over the background prediction, we set limits on black 
hole and string ball production rates. The following systematic uncertainties are taken into 
account in the limit setting procedure. The uncertainty in the background due to the choice 
of the fit function and its normalization rises with Sj from 1% to as much as 200% at very 
high values of St where the background extrapolation is unreliable, owing to insufficient data 
in the control regions. The uncertainties in the signal include the 8% uncertainty due to the 
jet energy scale, which is known to ~ 2% |2l) : a 6% uncertainty in the signal due to the PDF 
choice, as determined using the prescribed PDF4LHC recipe |37); and the 4.4% uncertainty 
in the integrated luminosity [38 1. As a result, the total systematic uncertainty in the signal is 
calculated to be 10%. 

For each set of model parameters, a test statistic S/ \/S + B, where S and B are the numbers 
of signal and background events, is used to choose an optimal combination of minimum Sj 
and multiplicity. Limits are then set using a modified frequentist CL S method Il39l l40l with a 
Poisson likelihood of the observed number of events, given the predicted background multi- 
plied by the likelihoods of a set of measurements of the nuisance parameters that are related to 
various systematic uncertainties, modelled by log-normal distributions. Counting experiments 
are performed to set a 95% confidence level (CL) cross section upper limit for each model used 
in this analysis. These limits can be interpreted in terms of lower mass limits on black holes 
(figure |4j and string balls (figure [5]) that range from 4.3 to 6.2 TeV. The mass limit plots show 
lower mass limits for a number of benchmark models as a function of the fundamental Planck 
scale, Mq. The areas below each curve are excluded by this analysis. 

For quantum black holes, which are characterized by a low final-state multiplicity N, the limits 
come from the N > 2 samples. As the N > 2 sample largely overlaps with the sample used for 
the background shape determination (N = 2), we use the N > 2 sample only to set limits on 
quantum black holes with masses above the range used for the background fit, as can be seen in 
figure]!] Note that the n = 1 case for quantum black holes corresponds to the RS black holes 0. 
In this case, Md is the Planck scale times the exponential factor coming from the warping of the 
anti-deSitter space, and is expected to be of the order of the electroweak symmetry breaking 
scale, similar to the fundamental Planck scale in the ADD model. The limits on the quantum 
black holes mass are shown if figure [5] All other benchmark model limits set in this paper 
correspond to the ADD model. 

To extend the scope of this search, model-independent limits on the cross section times the ac- 
ceptance (A) are computed for high-Sj inclusive final states for N > 3, 4, 5, 6, 7, 8, 9, and 10, as 
a function of minimum Sj (figures |6j[7|). The intersection of these limits with theoretical pre- 
dictions for the cross section within the fiducial and kinematic selections used in this analysis 
could be used to constrain other models of new physics resulting in energetic, multiparticle 
final states. These model-independent limits on the cross section times acceptance are as low 
as 0.2 fb at 95% CL for minimum Sj values above ~4.5 TeV, where no data events are observed. 

6 Summary 

A search for microscopic black holes and string balls at the LHC has been conducted, using 
a data sample corresponding to an integrated luminosity of 12.1 ± 0.5 fb -1 of y/s = 8 TeV pp 
collisions collected with the CMS detector at the LHC in 2012. Comparing the distributions 
of the scalar sum of the transverse momenta of all the final-state objects in data events with 
those from background, new model-independent limits are set that can be used to constrain a 
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Figure 4: The 95% CL lower limits on the semiclassical black hole mass derived from the upper 
95% CL limits on cross section times branching fraction as a function of the fundamental Planck 
scale Md, for various models. The areas below each curve are excluded by this search. Top 
left: BlackMax black hole models without the stable remnant. Top right and bottom row: 
Charybdis black hole models with or without the stable remnant. The number n of extra 
dimensions is labelled accordingly. 
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Figure 5: (Left) The cross section upper limits at 95% CL from the counting experiments opti- 
mized for various string ball parameter sets (solid lines) compared with predicted signal pro- 
duction cross section (dashed lines) as a function of minimum string ball mass. Here, Mp is 
the multidimensional Planck scale, Ms is the string scale, and g$ = 0.4 is the string coupling. 
(Right) Lower quantum black hole mass limits at 95% CL as functions of the fundamental 
Planck scale Md for various QBH black hole models with a number n of extra dimensions from 
one to six. 



wide variety of models. With this search, semiclassical and quantum black holes with masses 
below 4.3-6.2 TeV are excluded in the context of a number of benchmark models. Stringent 
limits on black hole precursors - string balls - are also set. These limits extend significantly the 
previously probed regime of black hole production at hadron colliders and represent the most 
restrictive exclusions on these objects to date. 
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Figure 6: Model-independent 95% CL cross section times acceptance (A) upper limits for count- 
ing experiments with Sj > S™ m as a function of S™ m for events with multiplicity: (Top left) 
N > 3, (top right) N > 4, (bottom left) N > 5, and (bottom right) N > 6. The blue solid 
(red dotted) lines correspond to an observed (expected) limit for nominal signal acceptance 
uncertainty of 10%. The green (dark) and yellow (light) bands represent one and two standard 
deviations from the expected limits. 
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Figure 7: Model-independent 95% CL cross section times acceptance (A) upper limits for count- 
ing experiments with Sj > S™ m as a function of S™ m for events with multiplicity: (Top left) 
N > 7, (top right) N > 8, (bottom left) N > 9, and (bottom right) N > 10. The blue solid 
(red dotted) lines correspond to an observed (expected) limit for nominal signal acceptance 
uncertainty of 10%. The green (dark) and yellow (light) bands represent one and two standard 
deviations from the expected limits. 
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